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ABSTRACT 

It has recently been discovered that some, if not all, classical novae emit GeV gamma-rays during 
outburst. Despite using an unreliable method to determine its distance, previous work showed that 
nova V1324 Sco was the most gamma-ray luminous of all gamma-ray-detected novae. We present 
here a different, more robust, method to determine the reddening and distance to V1324 Sco using 
high-resolution optical spectroscopy. Using two independent methods we derived a reddening of 
E(B — V) = 1.16 ± 0.12 and a distance limit of r D > 6.5 kpc. This distance is > 40% greater 
than previously estimated, meaning that V1324 Sco has an even higher gamma-ray luminosity than 
previously calculated. We also use periodic modulations in the brightness, interpreted as the orbital 
period, in conjunction with pre-outburst photometric limits to show that a main-sequence companion 
is strongly favored. 

Subject headings: novae, cataclysmic variables, stars: distances 


1. INTRODUCTION 

Classical nova are the result of a thermonuclear run¬ 
away taking place on the surface of a white dwarf and 
are fueled by matter accreted on to the white dwarf from 
a companion star. These outbursts give rise to an in¬ 
crease in luminosity and eject between ~ 10~ 3 —10~' Mq 
of material at velocities > 10 3 km s _1 (Prialnik 1986; 
Gallagher & Starrfield 1978; Shore 2012). 

Nova outbursts have also been detected in the GeV 
gamma-rays regime with Fermi Gamma Ray Space Tele¬ 
scope (see e.g. Cheung et al. 2010, 2012a,b; Hays et al. 
2013; Cheung et al. 2015). The likely source of 
these gamma-ray photons is strong shocks, which can 
generate relativistic particles via the diffusive shock 
acceleration mechanism and the relativistic particles 
can then generate gamma-rays (Metzger et al. 2014; 
Blandford & Ostriker 1978). But the details of shocks 
and shock acceleration in novae still remains a poorly un¬ 
derstood issue, despite the potential for insight into the 
nature of these high-energy events (Metzger et al. 2015). 

The first nova to be detected by Fermi was V407 
Cyg, and it received considerable attention (Abdo et al. 
2010; Chomiuk et al. 2012; Orlando & Drake 2012; 
Nelson et al. 2012). Given that V407 Cyg has a Mira gi¬ 
ant secondary with a strong wind (a member of the sym¬ 
biotic class of systems), a model to explain the gamma- 
rays was proposed wherein a shock was generated as the 
nova ejecta interacted with the dense ambient medium 
surrounding the red-giant companion. 

This model, however, failed to explain V1324 
Sco, V959 Mon, V339 Del (Ackermann et al. 2014), 
V1369 Cen (Cheung et al. 2013), and Nova Sgr 
2015b (Cheung et al. 2015), all of which are novae that 
were detected by Fermi but lack a detectable red- 
giant companion (hereafter referred to as non-symbiotic 
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gamma-ray-detected novae which, for the sake of brevity, 
we will refer to as NGDN). This non-detection of 
a red-giant companion implies that these novae have 
main-sequence companions with low-density circumstel- 
lar material. We note that, while it is possible for 
these novae to have high density circumstellar mate¬ 
rial despite not having a red-giant companion, no evi¬ 
dence has yet been found for dusty circumstellar mate¬ 
rial around cataclysmic variables (Harrison et al. 2013). 
The primary analysis for the NGDN was presented 
in Ackermann et al. 2014, which did not include V1369 
Cen or Nova Sgr 2015b. Therefore the discussion that fol¬ 
lows will only focus on the three non-symbiotic gamma- 
ray-detected nova discussed in Ackermann et al. 2014. 

While no direct explanation was proposed by 
Ackermann et al. 2014 to explain the gamma-ray pro¬ 
duction in V1324 Sco, V959 Mon, and V339 Del, they 
did suggest that these novae were similar enough to 
one another that they represented a homogeneous pop¬ 
ulation. Their evidence for a homogeneous population 
was that all three had similar gamma-ray spectra and 
light curves and were otherwise unremarkable novae. 
Ackermann et al. 2014 went on to suggest that it is only 
the close proximity of the NGDN that made them de¬ 
tectable by Fermi , and that other novae would be de¬ 
tected if only they were closer. However, as we show in 
this paper, the NGDN are not actually homogeneous, as 
one of them clearly stands out from the rest. 

VI324 Sco is located in one of the fields continu¬ 
ally observed by the Microlensing Observations in As¬ 
trophysics (MOA) group, and its outburst was initially 
detected in 2012 May as part of their high-cadence /- 
band photometry (Wagner et al. 2012). The initial de¬ 
tection showed a slow monotonic rise in brightness be¬ 
tween May 14 - May 31, followed by a very large (A I ~ 
6 mag) increase in brightness between June 1 - June 
3 (Wagner et al. 2012). Approximately two weeks after 
its optical brightening (2012 June 1, which we take to 
be the primary nova event) it was detected by the Fermi 
collaboration as a new transient source (Cheung et al. 
2012b), with gamma-ray emission lasting a further ~2 
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weeks (Ackermann et al. 2014). There are several things 
that make V1324 Sco unusual when compared to the 
other NGDN. 

To begin with, V1324 Sco is the most gamma- 
ray luminous, being at least ~ 2 times greater than 
the other gamma-ray-detected novae (including V407 
Cyg). The estimate for the gamma-ray luminosity, done 
by Ackermann et al. 2014, was based on distance esti¬ 
mations derived using the Maximum Magnitude Rate of 
Decline (MMRD) method, a technique that has been the 
subject of serious critiques in recent publications (see 
e.g. (Kasliwal et al. 2011; Cao et al. 2012). The MMRD 
derived distance used in Ackermann et al. 2014 was ~ 4.5 
kpc, derived using the relation of della Valle & Livio 
1995. It should be noted that Ackermann et al. 2014 
do emphasize the uncertainty associated with this value, 
explicitly stating that a possible dust event could com¬ 
promise their determination of the rate of decline in lumi¬ 
nosity, which would alter their derived distance. In this 
paper we find a distance limit 40% greater than that, 
meaning that the gamma-ray luminosity for VI324 Sco 
is even greater than initially calculated. 

Another unusual aspect of V1324 Sco is that it was 
never detected as an X-ray source. Along with dif¬ 
fuse shock acceleration, shocks also generate a tremen¬ 
dous amount of heat, raising the ejecta to ~ keV 
temperatures (Metzger et al. 2015, 2014). This hot 
shocked gas can be observed by looking for hard X- 
rays (Mukai & Ishida 2001; Mukai et al. 2008). A 
stronger shock should generate more heat, which im¬ 
plies that the gamma-ray luminosity and X-ray luminos¬ 
ity should correlate with one another. However, despite 
the extreme gamma-ray luminosity of V1324 Sco, it was 
never detected as an X-ray source (Metzger et al. 2014). 

V1324 Sco clearly stands out among the NGDN, and 
this paper is focused on placing better constraints on 
this unique nova by finding values for the reddening and 
a distance limit. We present the results from our analy¬ 
sis of both spectroscopic and near-IR photometric data 
of V1324 Sco to find a distance limit, reddening, and 
constrain the luminosity class of the companion star. In 
Section 2 we present the details of the data observations 
and reduction. In Section 3 we detail the methods used 
to derive a value for the reddening, and in Section 4 we 
provide similar details for our distance limit. In Section 5 
we discuss the constraints we can place on the progenitor 
system and in Section 6 we conclude by deriving a new 
lower limit on the gamma-ray luminosity using our new 
distance limit. 

2. OBSERVATIONS 

To determine the reddening we analyzed two high- 
resolution spectroscopic observations of V1324 Sco; a 
very early time Very Large Telescope (VLT) observa¬ 
tion, and a late time observation using the Magel¬ 
lan telescope. The VLT observation was taken from 
the archive 4 , with the original observations reported by 
Wagner et al. (2012). The Magellan observations were 
obtained by the authors. 

2.1. VLT UVES 

4 Based on observations made with ESO Telescopes at the 
Paranal Observatory under programme ID 089.B-0047, PI: Feltzing 


Table 1 

Best Fit Reddening Values For DIBs 


DIB A (A) 

EW (mA) 

E(B - V ) 

5705.1 

99 ± 9 

1.01 ±0.11 

5780.5 

581 ± 40 

1.14 ±0.08 

5797.1 

155 ± 16 

0.86 ±0.09 

6196.0 

68 ± 2 

1.39 ±0.04 

6613.6 

200 ± 8 

0.94 ±0.04 


Spectra were taken on 2012 June 4.1 UT (+3.1 days 
after main outburst) by Wagner et al. 2012 using the Ul¬ 
traviolet and Visible Echclle Spectrograph (UVES) in¬ 
strument on the VLT (Dckker et al. 2000). The spec¬ 
tra were taken as part of the follow-up conducted by 
the MOA group. Copies of the data were made publicly 
available on the ESO archive and were obtained by the 
authors. 

Observations were made in dichroic mode, with the 
blue arm centered at 4370 A (spanning 3600 — 4800 A) 
and the red arm centered at 7600 A (covering 5600 — 9300 
A), taken using two CCDs that have a chip gap between 
7550 — 7650A. The slit width was 1.00" with an average 
seeing of 0.8", giving a resolution of R ~ 40, 000 for the 
blue arm and 40,000 — 50,000 for the red arm. Four 
1800 second integrations were taken with lxl binning, 
with each integration having an average signal-to-noise of 
S/N « 30 per pixel for the blue arm and S/N ss 110—175 
per pixel for the red. 

The reduction was undertaken using the standard ESO 
Reflex data reduction pipeline, which includes Hatfield 
correction, bias subtraction, cosmic-ray removal, spectral 
extraction, and wavelength calibration using comparison 
spectra of a ThAr lamp (see Sacco et al. 2014 for details 
on the data reduction procedure). The four individual 
spectrum frames were combined using the IRAF routine 
ndcombine , giving a final S/N per pixel of ~ 60 in the 
blue and S/N per pixel of « 220 — 350 in the red. 

2.2. Magellan MIKE 

Further observations were made on 2012 July 16.1 UT 
(+45.1 days after main outburst) using the MIKE instru¬ 
ment (Bernstein et al. 2003) on the 6.5 meter Magellan 
Clay telescope. 

The MIKE instrument also has two arms, blue (3350 — 
5000 A) and red (4900—9500 A). Two 300 second integra¬ 
tions were taken with a 0.7" slit, giving a final resolution 
of R ss 40,000 in the blue and R ~ 30,000 in the red. 
CCD binning was 2x2 for all integrations, giving an 
average S/N per pixel of ~ 30 in the blue and S/N per 
pixel of ~ 170 in the red. 

The reduction for the MIKE spectra was done using 
the Carnegie Python tools (CarPy) 5 , which provides a 
simple pipeline data reduction procedure. Milky flats, 
Quartz flats, Twilight flats,ThAr comparison lamps and 
bias frames were all taken during the observation run and 
utilized in the data reduction pipeline. 


5 http://code.obs.carnegiescience.edu/mike 
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Figure 1. Five DIBs that were used to determine reddening. The largest contributor of uncertainty in the EW of the DIB features was 
determination of the continuum flux. This is why, in spite of the poor S/N, the 5705.1 A feature has a comparable uncertainty to the other 
features. 



Figure 2. The best fit E(B — V) values for all five DIBs that 
were analyzed (see Table 1 for the specific values). The redden¬ 
ing value for each DIB was determined using the best-fit val¬ 
ues of Friedman et al. (2011). The average reddening value is 
E(B — V) = 1.07 (the solid horizontal blue line) with an uncer¬ 
tainty of ±0.20 (the dashed horizontal red lines). The green band 
corresponds to the reddening value, E{B — V) = 1.23 ± 0.15, de¬ 
rived using the EW of Na and K absorption features. The width 
of the band corresponds to the uncertainty in the derived value. 

3. REDDENING 

3.1. Reddening Measurement from DIBs 

We use the Equivalent Width (EW) of Diffuse In¬ 
terstellar Bands (DIBs) to measure the reddening, a 
method that has already been used in the context of 
novae (see e.g., Shore et al. 2011; Munari et al. 2012). 
Friedman et al. (2011) found an empirical relationship 
between the EW of eight strong optical DIBs and the 


reddening along a given sight line. In V1324 Sco, one of 
these eight DIB features (5487.7 A) had very low S/N 
(in both VLT and MIKE spectra) and was not used, and 
two features (6204.5 A and 6283.8 A) were discarded due 
to overlapping telluric features. The remaining five DIBs 
were used to find a value for E(B — V). All of the EW 
values were measured using the VLT spectrum (see Fig¬ 
ure 1) as it had a much higher S/N in the red compared 
to the MIKE spectrum. The specific DIB EW values are 
given in Table 1, and the derived E(B — V) as a function 
of wavelength is plotted in Figure 2. All EW uncertain¬ 
ties were determined by iterative fitting solutions using 
the IRAF tool splot to determine the spread of potential 
EW values. The primary uncertainty in determining the 
EW values was setting the continuum flux level. From 
our derived values for the reddening we took an error 
weighted average and found E{B — V) = 1.07 with a 
la uncertainty on the error weighted average of ±0.20. 
The uncertainty for the individual E(B — V) values de¬ 
rived for each DIB included both the uncertainty in the 
fit parameters from Friedman et al. (2011) as well as the 
standard deviation in measured EW values. The final 
uncertainty in the error weighted average was dominated 
by the spread in derived E(B — V ) values for different 
DIB features. 

3.2. Reddening Measurement from Na and K 
Absorption Features 

We used a second, independent, method to determine 
the reddening. This method utilized the empirical re¬ 
lationship found by Munari & Zwitter (1997), which re¬ 
lates reddening to the total equivalent width of the Na I 
D absorption lines (at 5889.9 A and 5895.9A) and K I ab¬ 
sorption line (at 7698.9 A). The material along the line of 
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Figure 3. Sodium D (5889.95 A) and potassium (7698.96 A) interstellar absorption features that were used to derive a value for the 
reddening. The specific absorption components from Table 2 are marked as black dashed lines. We plot both the UVES and Magellan 
spectra for the sodium to show that these features are unchanging in time, meaning they are interstellar and not associated with the nova. 
Note that we only use the potassium absorption feature to calculate reddening when the sodium feature is saturated, which is why we are 
only showing the potassium in the region where the sodium is saturated. 


sight to an object will have some characteristic features 
due to interstellar absorption, and Munari & Zwitter 
(1997) calibrated two of these features to determine the 
amount of reddening each feature contributes. The sum 
of reddening from individual absorbing features gives the 
total reddening. 

Munari & Zwitter (1997) found that the Na I D fea¬ 
tures are ideal for tracing reddening at low column den¬ 
sities, but saturates at high column densities. In these 
instances the K I feature, which does not saturate as eas¬ 
ily, can be used to determine the reddening. In our case, 
we only needed to utilize the K I for the two Na I D 
features that were saturated (at —7 and +3 km/s). 

To insure that the absorption features were the result 
of interstellar clouds, and not the nova itself, we com¬ 
pared the UVES and MIKE spectra. If the features were 
from the nova, we would expect them to change over 
time; as can be seen in the Na I D plot, which shows both 
the UVES and MIKE spectra, the features remained con¬ 
stant. We did not plot the MIKE spectra for the K I 
features, as the S/N was very low. 

To avoid potential contamination of the lines with tel¬ 
luric features we used an archived telluric divider to re¬ 
move telluric features from the UVES spectrum. The 
EW, derived reddening, and velocity of the absorption 
features that we used to find the total reddening are given 
in Table 2 and Figure 3. Although the uncertainty re¬ 
sulting from sum of the measured quantities is ±0.09, 
the intrinsic scatter that Munari & Zwitter (1997) found 
in their calibration was 0.15 for large reddening values 
( E(B — V) > 0.4). As a result, we take 0.15 to be the 
uncertainty. 

The reddening value found using this technique, E(B — 
V) = 1.23 ± 0.15, is consistent (within margin of uncer¬ 
tainty) with the reddening derived using the DIB fea¬ 
tures. 

We combine these two independent reddening mea¬ 
surements by taking an error weighted average of the 
two, which yields our final reddening value of E(B — 
V) =1.16 ±0.12. 

4. DISTANCE 


Table 2 

Reddening Values from Na I D & K I Absorption Lines 


Absorption 
Wavelength (A) 

EW (mA) 

to 

1 

LSR Radial 
Velocity (km s — - 1 ) 

5886.94 

246 ± 12 

0.08 ± 0.01 

-136 

5887.15 

127 ± 6 

0.04 ± 0.01 

-126 

5888.06 

148 ± 7 

0.05 ± 0.01 

-79 

5888.42 

234 ± 12 

0.08 ± 0.01 

-61 

5888.82 

171± 9 

0.05 ± 0.01 

-41 

7698.37 

103 ± 5 

0.39 ± 0.02 

-7 

7698.62 

140 ± 7 

0.54 ± 0.02 

3 

Total 


1.23 ± 0.09 



In order to derive distance we used the reddening value 
found in Section 3 in conjunction with a 3D Galactic cen¬ 
ter reddening map from Schultheis et al. (2014), found 
using data from the Vista Variables in the Via Lactae 
(VVV ) survey. The VVV survey is a ESO large pro¬ 
gram using the 4-meter VISTA telescope to take near-IR 
photometry (0.9 — 2.5^m) of 520 square degrees towards 
the Milky Way Bulge to characterize variable sources. 
The reddening map was one of the byproducts of the 
exquisite VVV photometric dataset. The 3D map gives 
E{J — K) reddening values as a function of: Galactic 
longitude (in intervals of 0.1°), Galactic latitude (in in¬ 
tervals of 0.1°), and radial distance (in intervals of 0.5 
kpc, extending out to 10 kpc). 

To make use of the map we needed to transform our 
E(B — V) value into an E(J — K) value given in the 
reddening map. To do this we found a coefficient, 7 , 
such that E(J — K) = 7 E(B — V). The value of 7 was 
determined using the VVV survey color transforms from 
Table 1 in Saito et al. (2012). These color transforms 
give the relative extinction for the VVV filter system in 
terms of E(B — V), assuming a standard extinction law 
of Cardelli et al. (1989). From these correction values we 
found a 7 value of 0.502, which gives E{J — I\) = 0.58 ± 
0.06 for our reddening value of E(B — V) =1.16 ± 0.12. 

We used the average of the four points in the redden¬ 
ing map closest to the coordinates of V1324 Sco (RA = 
17:50:53.90, Dec = —32:37:20.5) for our analysis. Fig¬ 
ure 4 shows E(J — K ) versus radial distance td for the 
reddening map values, as well as our derived value for 
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Figure 4. The green data points give the (average) E(J — K) 
reddening values as a function of distance at the position of V1324 
Sco, taken from the 3D reddening map of Schultheis et al. (2014). 
The red diamond indicates the derived reddening. The blue dashed 
lines show the lcr extent of the derived reddening; it is clear that 
this value is consistent with all reddening values past 6.5 kpc. 

V1324 Sco’s E(J — K) reddening. 

The dashed blue line in Figure 4 shows the lcr max¬ 
imum for our derived E(J — K) reddening (0.62); it is 
clear that this value is consistent with any distance >6.5 
kpc. Because of this degeneracy in E{ J — K) reddening 
we can only place a lower limit of 6.5 kpc on the distance 
to V1324 Sco. 


5. PROGENITOR 

If VI324 Sco had a giant companion like V407 Cyg 
then the wind from the companion could explain its ex¬ 
ceptionally high gamma-ray luminosity. To determine if 
this is a feasible explanation we need to determine if the 
progenitor system to V1324 Sco has a giant secondary 
star. The strongest photometric constraints on the sec¬ 
ondary star comes from the VVV Survey (Minniti et al. 
2010), which obtained several epochs of near-IR photom¬ 
etry on the target field between July 2010 and September 
2011. We searched the catalog for all stars with non-zero 
AT-band magnitude in a 15’xl5’ cutout centered on the 
coordinates of V1324 Sco (N « 25,000). To qualify as a 
star an object had to have a pStar value > 0.9 (see VVV 
documentation for a description of the pStar variable). 
From this sample we found that 99% of the sources had 
rriK < 16.626; given that V1324 Sco was not detected by 
the VVV survey, we use this value as a limiting magni¬ 
tude. 

Table 1 in Saito et al. (2012) gives the filter specific ex¬ 
tinction for the VVV K band as Ak/E(B — V) = 0.364, 
and we can use our derived reddening value to obtain a 
AT-band extinction of ~ 0.42 mags. Without a proper 
upper bound on the distance to VI324 Sco we cannot 
place an upper limit on the distance modulus. How¬ 
ever we can say that, if V1324 Sco were at a distance 
of > 9.5 kpc, then the distance modulus would be 
ss 15.3, meaning that the absolute AT-band limit would 
be bright enough for the companion to be a very faint 
giant (M-k > 1-3, corresponding to a spectral type G5 
giant (Covey et al. 2007)). However, if V1324 Sco is at 
the galactic center (~ 8.5 kpc), then it must have a dwarf 


o 

y —I 

I 


I 


M K 

O 


in 


1 10 100 1000 10 4 10 5 10 6 
Orbital period (hours) 

Figure 5. Plot showing the distribution of novae as a function 
of absolute X-band magnitude and binary orbital period, taken 
from Darnley et al. (2012). V1324 Sco is marked as the purple 
transparent line, indicating that we know the period ~ 1.6 hours, 
but we cannot place a limit on the absolute K -band magnitude. 
The width of the purple line corresponds to the fact that the period 
could either be 1.6 or 3.2 hours. The colors on the plot correspond 
to the luminosity class of the secondary: blue is main sequence, 
green is sub-giant, and red is giant. 

companion. 

While this limit isn’t stringent enough to rule out a 
giant companion, by analyzing the period of VI324 Sco 
we show that the a main sequence companion is strongly 
favored. 



5.1. Constraints from the Orbital Period 

Darnley et al. (2012) created a classification system for 
novae according to their orbital period which is a proxy 
for luminosity class of the secondary star. We utilize this 
classification scheme to help constrain the companion to 
V1324 Sco. 

Measurements of the binary orbital period were found 
in the original photometry from the MOA group, who de¬ 
tected periodic modulations in the brightness of V1324 
Sco (Wagner et al. 2012). These modulations were on 
the order of ~ 0.1 mag, with a period of ~ 1.6 hours. 
However, it is possible that there cllipsodial variations 
in the photometric light curve, which would give a sec¬ 
ondary minimum and maximum. As a result, we consider 
both 1.6 and 3.2 hours for the period in our analysis. 

Interestingly, if the period is 1.6 hours, V1324 Sco 
would be below the period gap, meaning that the an¬ 
gular momentum loss is driven by gravitational radia¬ 
tion Knigge 2011. Only a handful of novae fall below the 
period gap (Wagner et al. 2012), and it has been pos¬ 
tulated that these systems may have a different type of 
mass transfer mechanism that takes place in this period 
range (Uthas et al. 2010). 

Using the value for the period within the framework 
of Darnley et al. (Figure 5) we see that V1324 Sco (the 
purple dashed line) most likely falls into the region of blue 
points, indicting that it has a main sequence companion. 
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6. CONCLUSIONS 

Here we present results of the analysis of optical spec¬ 
troscopy and near-IR photometry of the nova V1324 Sco 
to determine the reddening and a limit on the distance. 
We derived a reddening value using two independent 
methods and found E(B — V) = 1.16 ± 0.12. Using this 
reddening value we found a distance lower limit of limit 
of ru > 6.5 kpc to V1324 Sco, which stands in contrast 
to the values of Ackermann et al. (2014), determined us¬ 
ing the unreliable MMRD technique, of 4.5 kpc. We 
also found that, while we cannot rule out a giant com¬ 
panion, the short period of V1324 Sco strongly favors a 
main sequence companion. V1324 Sco was already the 
most gamma-ray luminous nova when using the values 
from Ackermann et al. (2014); using our new distance 
limit of rjj > 6.5 kpc, we find that the gamma-ray lu¬ 
minosity is at least a factor of ~ 2.1 times greater than 
initially calculated by Ackermann et al. (2014). We en¬ 
courage deep follow-up photometry of V1324 Sco in order 
to place greater constraints on the companion. 
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